Interplay between magnetic and lattice excitations in BiFeOa 
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We have performed Raman scattering investigations on the high energy magnetic excitations in a 
BiFeOa single crystal as a function of both temperature and laser excitation energy. A strong feature 
observed at 1250 cm -1 in the Raman spectra has been previously assigned to two phonon overtone. 
We show here that its unusual frequency shift with the excitation energy and its asymmetric tem- 
perature dependent Fano lineshape reveal a strong coupling to magnetic excitations. In the same 
energy range, we have also identified the two-magnon excitation with a temperature dependence 
very similar to a-Fe203 hematite. 
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Recently, multiferroics have attracted much attention 
because they exhibit a strong coupling between magnetic 
and ferroelectric orders X They are good candidates for 
controlling the magnetic properties by electric field and 
vice versa£ making these materials promising for poten- 
tial applications in spintronics and information storage. 
BiFeOa (BFO) is one of the most studied multiferroic due 
to both its strong ferroelectric polarization 3 -^ and anti- 
ferromagnetic order at room temporal lire'' ~ s . The ori- 
gin of the magnetoelectric effect in BFO remains an open 
issue. In order to elucidate this question, investigating 
the interplay between lattice vibrations and magnetic ex- 
citations appears very useful because lattice distortions 
strongly affect the ferroelectric polarization and therefore 
its coupling to magnetic order. 

THz spectroscopies have recently shown to be very 
efficient tools for probing the coupling between the 
lattice and magnetic excitations. Raman 9 - ' 10 ' 11 and 
infrared measurement a 12 ' 13 combined with theoretical 
approaches^ have permited to identify most of the vi- 
brationnal normal modes in BFO. At the same time 
low-frequency magneto-optical resonances in the dielec- 
tric susceptibility of multiferroic compounds have been 
observe d 15 ' 16 ] 17 ] 18 and recently the cycloidal magnetic 
modes have been optically detected in BFO 19 ' 20 . How- 
ever, probing magnetic excitations and in particular their 
coupling to the lattice is still a challenge in BFO because 
we are still lacking experimental evidences for the entan- 
glement with magnetic excitation through the phonon 
behaviors. 

Here, we present Raman measurements on BFO single 
crystals aiming to address the interplay between lattice 
and magnetic orders. Among the two-phonon overtone 
features previously detected at high energy (above 700 
cm-1), we show that the mode located at 1250 cm" 1 
is much more complex than a simple double vibrational 
mode and exhibits a strong magnetic character. Indeed 
we reveal that the 1250 cm" 1 mode depends on the laser 
excitation energy and exhibits an asymmetric lineshape 
controlled by the spin dynamic of BFO. We have also 
detected in the same energy range, the signature of a 



two magnon mode near 1050 cm" 1 , its temperature- 
dependence is in agreement with the one observed in a- 
Fe2 03 hematite . 

BiFeOa single crystals were grown in air using a Bi2 03- 
Fe2 03 flux technique^ They are typically of millimeter 
size. Below a high Curie temperature (T c ^1100 K) BFO 
has a rhombohedrally distorted perovskite structure with 
a spontaneous electrical polarization (between 50 and 
100 /iC/cm 2 in magnitude) along the [111] direction in 
the pseudo-cubic configuration^. Below the Neel temper- 
ature (T/v~640 K)2il, BFO exhibits a G-type antiferro- 
magnetic order subjected to long range modulation as- 
sociated with a cycloidal spiral with a period length of 
62 nm 8 -. The spiral propagation direction is along [10-1] 
with a spin rotation within (-12-1) plane. 

Raman measurements have been performed using sev- 
eral laser excitation lines between 1,92 to 2.41 eV from 
a mixed Ar + -Kr + gas laser. The BFO crystals were 
mounted in vacuum (10 _6 mbar) on the cold finger of 
a liquid helium cryostat. The Raman spectra were 
recorded between 7 and 300 K using a triple grating spec- 
trometer (JY-T64000) equipped with a nitrogen cooled 
CCD detector. The spectrometer was in substractive 
configuration with a resolution of 0.5 cm" 1 . The Raman 
spectra have been obtained using unpolarized excitation 
light. 

The assignment of the vibrationnal modes under 
700 cm" 1 is well known] 9 ' ' n Group theory predicts 13 
Raman active phonons for the R3C symmetry of BFO: 
T = 4Ai + 9E. Above 700 cm -1 , one should only detect 
the two-phonon Raman peaks. Two-phonon Raman scat- 
tering is a second order process involving two phonons. 
In contrast to one-phonon scattering process there is no 
restriction on the phonon wavevector but the sum of the 
two phonon wavectors has to be zero. If the two phonons 
are identical, the two-phonon Raman peak is called an 
overtone. 

Figure Q] shows the Raman spectra in the 1000- 
1300 cm" 1 range as a function of the laser excitation 
energy. We focus here on the three main peaks at 1085, 
1140 and 1250 cm" 1 (indicated by arrows in Fig.[l|. The 
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FIG. 1: Raman spectra of BiFe03 single crystal in the 900- 
1500 cm -1 frequency range at 20K. The spectra labeled A to 
H were obtained from 2.71, 2.65, 2.6, 2.54, 2.41, 2.34, 2.18, 
1.92 eV excitation laser lines, respectively. The Raman peaks 
of are indicated by arrows. Dashed lines are guides for the 
eyes. All spectra have been corrected for the spectral response 
of the spectrometer. 



Raman peak intensities are strongly enhanced under the 
blue (2.34 eV) excitation line. These enhancements or 
resonance come from the matching of the incident pho- 
ton energy with the absorption edge of BFO i 10 i 21 As the 
incident photon energy is close to the absorption edge, 
the probability of creating electron-hole pairs diverges 
and the Raman scattering efficiency is greatly enhanced. 

The three peaks have been assigned to a two phonon 
Raman scattering process^ In particular, the promi- 
nent peak at 1250 cm -1 has been previously reported 
by Fukumura et al.. Its sudden disappearence has been 
observed around 600 K-700 K close to pointing out 
it's probably magnetic origin^ More recently, Ramirez 
et al. have assigned this peak to a pure two phonon pro- 
cess (2x[E g (TO)= 620 cm^ 1 ])^ However, the decrease 
of its integrated intensity up to TV and its constant be- 
haviour beyond Tjv advocate in favor of a strong coupling 
between the ferroelectric and magnetic subsystems. 

In Fig.[TJ the two phonon modes locations at 1085 and 
1140 cm -1 remains constant in frequency as the excita- 
tion energy is tuned. On the contrary, the frequency of 
the 1250 cm -1 peak increases when the excitation energy 
is tuned from 2.71 to 2.41 eV and then decreases. The 
frequency shift is about 20 cm -1 as shown in Fig. HJb). 
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FIG. 2: (a) Intensities normalized to its maximum value and 
(b) wavenumber shifts as a function of the laser energy for the 
1085 (diamond), 1140 (star) and 1250 cm -1 (square) peaks. 
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FIG. 3: Temperature dependance of the two-magnon and 
1250 cm -1 peaks obtained with the 2.34 eV laser line. Dashed 
lines are guides for the eyes. The peak at 1250 cm -1 has been 
fitted by the Fano equation (see text). The inset shows the 
temperature evolution of the asymmetric parameter q corre- 
sponding to the function normalized to its maximum value. 



To explain the different behaviour with the excitation en- 
ergy of the 1250 cm -1 overtone compare to the 1085 and 
1140 cm -1 one different scenarii can be put forward: 

(i) The peak at 1250 cm -1 is broad and might be con- 
stituted of two peaks with two different resonant behav- 
iors with the excitation energy thus explaining the ob- 
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FIG. 4: Temperature dependence of the normalized two- 
magnon frequencies for BiFe03 (square) and hematite a- 
Fe2 03 (diamond). 



served shift of the 1250 cm -1 peak in Fig. [TJ The fre- 
quency of these two peaks located at 1228 and 1288 cm -1 
can be deduced from a gaussian fit deconvolution of the 
broad 1250 cm" 1 peak. The 1228 and 1288 cm" 1 peaks 
should result respectively from an overtone of a peak 
around 615 and 645 cm -1 . However, only one peak has 
been detected around 620 cm -1 and we are led to con- 
clude that the peak observed at 1250 cm -1 is unlikely 
composed of two peaks. 

(ii) The Raman band is associated with a lumines- 
cence process assisted by magnons or phonons. Strong 
exciton-phonon interactions lead to new exciton-phonon 
bound states which causes phonon-assisted exciton tran- 
sitions. For luminescence, the transitions via the new 
bound states induce a continuous background which in- 
terferes with the usual discrete transitions imply in the 
luminescence process. This would be the origin of Fano 
lincshape in the phonon-assisted luminescence process. 
Since luminescence is produced by thermalized electron- 
hole pairs, the emmitted photons have no correlation 
with the excitation lines. As a result, the frequency of the 
luminescence peak is independant of the excitation en- 
ergy and increases monotonically with the incident pho- 
ton energy. The frequency maximum of the 1250 cm -1 
peak in Fig. (Hb) is not consistent with this interpreta- 
tion. 

(iii) Such a behaviour is due to the coupling of the 
peak at 1250 cm -1 with the antiferromagnetic order as 
suggested by Fukumura et al. and Ramirez et ai. 

In order to analyse and identify the different contri- 
butions of the three main peaks detected in the 1000- 
1300 cm -1 range of the Raman spectra, we have pre- 
formed Raman measurements as a function of tempera- 



ture (from 20 to 300 K) as shown in Fig. [3l The three 
main peaks exhibit a small softening in energy (about 
5 cm -1 ), gradually broaden and decrease in intensity 
as the temperature is raised. The peak at 1250 cm -1 
presents an asymmetric temperature dependent linc- 
shape characteristic of a strong coupling between the two 
phonon mode and the spin excitation spectrum as we will 
see just below. 

In order to quantify the Fano lineshape we have fit- 
ted the Raman peak located at 1250 cm -1 by the Fano 

equation^ I{uj) = I ^ g+ (u ,_^ o) ^ where ujq is the peak 

position without interaction, q is the Fano asymmetric 
parameter and T is the width at half maximum. The q 
parameter serves as a measure of the coupling between 
the two phonon state and a continuum of states that we 
have to identify. 

The evolution of the normalized asymmetric parameter 
q as a function of the temperature is plotted in the inset 
of Fig. [3] A first glance on the temperature evolution 
of the q parameter reveals two distinct regimes: q(T) 
exhibits a sharp step close to T — 150 K. 

Recent measurements on spin waves in BFO 
have shown a spin-reorientation close to these two 
temperatures! 19 ! 20 The q parameter is thus sensitive to 
the spin-rcoricntation which points out the intimate re- 
lationship between the double phonon peak at 1250 cm -1 
and the spin sublattice excitations. From 10 K to 150 K 
the asymmetric parameter q remains constant and keeps 
its high value indicating that the two-phonon mode is 
strongly coupled to the spin excitations in this temper- 
ature range. Around 150 K, the spin-phonon coupling 
is modified by the reorientation of the spins and sharply 
decreases as shown by the decrease of the q parameter. 
An additional feature acures at 200 K where another 
spin-reorientation has been observed^ - Above 200 K the 
asymmetry parameter decreases continuously as a mir- 
ror of the softning of the spin-phonon coupling due to 
thermal excitations. 

In Fig. [21 the two phonons peaks at 1085 and 
1140 cm -1 disappear at 140 K leaving a broad and weak 
peak. From 140 K to 300 K this broad peak shifts from 
1110 cm -1 down to 1050 cm -1 . The large energy soften- 
ing (about 60 cm -1 ) of the peak with temperature is un- 
usual compared to the one of the phonons (about 5 cm -1 
as mentionned before). We suspect that this large fea- 
ture is related to a magnetic excitation and thus sensitive 
to Tjv. 

In order to establish the origin of this broad peak, 
we have studied the temperature dependence of its en- 
ergy. In Fig. [4] is plotted the normalized frequency shift 
versus the temperature normalized to T^. The broad 
peak contibution in Fig. [3] has been obtained substract- 
ing the gaussian fit contribution of the phonon peaks. 
We observe softening of the freqeuncy of the presumed 
two-magnon peak as a function of temperature quite 
consistent with the temperature dependance of the two- 
magnon peak previously measured in hematite a-Fe2C>3, 
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FeBC>3 and MnFz^i. For clarity, only the behaviour of 
the two-magnon peak in hematite has been plotted in 
Fig. [H The good agreement between the two curves 
strongly supports the assignment of the broad peak cen- 
tered at 1050 cm -1 (300 K) to a two-magnon peak. 

The two-magnon scattering involves two-magnons with 
equal and opposite wave vectors due to the momentum 
conservation similar to the two-phonon Raman scattering 
process The peak position related to the two-magnon 
scattering is expected equal to twice the zone boundary 
magnon energy and can be compared to the one measured 
in orthoferrites (RFe03 with R a rare-earth atom). In 
YFeC>3 and ErFeOs the two-magnon band is measured 
around 900 cm" 1 ^! 

We point out that the two-magnon signature should 
also be observed in infrared spectroscopy. In some crys- 
tal with a center of inversion the two-magnon excita- 



tion through the exchange dipole moment is not allowed 
but can be observed using the oscillation strength of the 
phonon. In BFO, no inversion symmetry is found be- 
tween the Fe 3+ ions which allow them to make significant 
contributions to the direct optical excitation of the two 
magnons. 

In summary, we have tracked the magnetic excitations 
of the high energy Raman spectra. The complex nature 
of the 1250 cm -1 feature have been revealed. This peak 
is an overtone of two phonon band with a frequency and 
a temperature dependance driven by its interaction with 
the magnetic subsystem. In addition, we have measured 
a weak and broad peak close to 1050 cm -1 whose 
temperature dependence is consistent with two-magnon 
scattering process. 
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